Hidden Fermi Liquid, Scattering Rate Saturation and Nernst Effect: 

Perspective 



aDMFT 



m 

o 

(N 

Oh 
< 

oo 



a 

o 



> 

00 

^* 

o 



X 



Wenhu Xu,^ Kristjan Haule/ and Gabriel Kotliar^ 

^Department of Physics and Astronomy, Rutgers University, 136 Frelinghuysen Rd., NJ 08854 

(Dated: April 30, 2013) 

We investigate the transport properties of a correlated metal within dynamical mean field theory. 
Canonical Fermi liquid behavior emerges only below a very low temperature scale Tfl - Surprisingly 
the quasiparticle scattering rate follows a quadratic temperature dependence up to much higher 
temperatures and crosses over to saturated behavior around a temperature scale Tsat- We identify 
these quasiparticles as constituents of the hidden Fermi liquid. The non-Fermi liquid transport 
above Tfl, in particular the linear-in-T resistivity, is shown to be a result of a strongly temperature 
dependent band dispersion. We derive simple expressions for resistivity. Hall angle, thermoelectric 
power and Nernst coefficient in terms of a temperature dependent renormalized band structure and 
the quasiparticle scattering rate. We discuss possible tests of the DMFT picture of transport using 
ac measurements. 

PACS numbers: 71.10.Ay, 71.27.+a, 72.10.Bg 



Fermi liquids [1] are good conductors. Quasiparticles 
{qp) with a mean free path much longer than their wave- 
length are responsible for the electric transport, and the 
resistivity vanishes quadratically at low temperatures. 
Landau theory is very robust and when reformulated in 
terms of a transport kinetic equation, it can be used to 
describe situations where Landau qps are strictly speak- 
ing not well defined, namely when the qp scattering 
rate is comparable to their energy, such as the electron- 
phonon coupled system above the Debye temperature [2] . 

The metallic state of many strongly correlated mate- 
rials is not described by Landau theory in a wide range 
of temperatures. Quadratic temperature dependence of 
resistivity occurs in a very narrow or vanishing range 
of temperatures. The interpretation of the resistivity in 
terms of the standard model of transport which is based 
on qps is problematic since it leads to mean free paths 
shorter than the {qp) DeBroglie wavelength as stressed 
by Emery and Kivelson [3]. The transport properties of 
these "bad metals" thus requires a novel framework for 
their theoretical interpretation. 

Dynamical mean field theory (DMFT) [4], provides 
a non-perturbative framework for the description of 
strongly correlated materials. It links observable quan- 
tities to a simpler, but still interacting, reference system 
(a quantum impurity in a self-consistent medium) rather 
than to a free electron system, hence it gives access to 
physical regimes outside the scope of Landau theory. 

In a broad temperature range, the single-site DMFT 
description of the one band Hubbard model at large U 
and finite doping, results in transport and optical prop- 
erties with anomalous temperature dependence [5-12], 
reminiscent of those observed in bad metals. Corre- 
sponding studies of half filled metallic systems [13-16] 
also reveal bad metallic behavior in a narrower temper- 
ature region since at high temperatures the resistivity is 
insulating like. 

Landau qps only emerge below an extremely low tem- 
perature, TpL, which is much lower than the renormal- 



ized kinetic energy or Brinkman-Rice scale Tbr ~ SW 
with S the doping level and W the bare bandwidth. Tbr 
is the natural scale for the variation of physical quanti- 
ties with doping at zero temperature [7, 17]. A recent 
comprehensive DMFT study of the Hubbard model with 
a semicircular bare density of states found that the trans- 
port properties above T^l are described in terms of re- 
silient qps with a strong particle-hole asymmetry [11]. 
This asymmetry arises from the asymmetric pole struc- 
ture in the self energy characterizing the proximity to the 
Mott insulator [18]. 

In this Letter we investigate the problem of bad metal 
transport. By expressing the DMFT transport coeffi- 
cients in terms of qp quantities we find several surpris- 
ing results: a) the qp scattering rate has a quadratic 
behavior for temperature much larger than Tfl and 
crosses over to a saturated behavior around Tgat- b) The 
temperature dependence of the transport coefficients is 
anomalous (in the sense that it does not reflect the T 
dependence of the qp scattering rate) and arises from 
the temperature dependent changes of the qp disper- 
sion near the Fermi level, c) The temperature depen- 
dence of the qp dispersion affects differently the diago- 
nal and off-diagonal charge and thermal transport co- 
efficients but the Mott relation [19, 20] is valid when 
Tfl<T< Tsat/2. 

We study the one-band Hubbard Hamiltonian on the 
two-dimensional square lattice with nearest neighbor 
hopping. 
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We set the full bare bandwidth W = 8t to W = 1 as the 
unit of energy and temperature, and present results for 
U/W = 1.75, for which the system is a Mott insulator 
at half-filling. The doping level of the metallic state is 
fixed dit S = 15% (n = 0.85). We use continuous time 
quantum Monte Carlo method (CTQMC) [21] and the 
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FIG. 1: Spectral function Ak(cj) along T-X-M-T in 
Brillioun zone at (a) T = 0.0075, (b) T = 0.025, (c) T = 0.1 
and (d) T = 0.125. (e) Local density of states, (f) Roots of 
Eq. 3 at different temperatures. 



The high temperature band breaks into two parts, which 
together with the upper Hubbard band forms the char- 
acteristic DMFT three-peak structure of local density of 
states (LDOS, Fig. 1(e)). The breakup of these bands 
also leads to the separation of the optical spectrum into 
a Drude peak and a mid-infrared feature, characteristic 
of many correlated systems, which provided the earliest 
experimental tests of the DMFT picture of correlated 
materials [5, 6, 13, 14] . 

There is always a dispersive qp feature in a ~ /c^T en- 
ergy window of the Fermi energy, uj^ denotes the root of 
Eq. 3 closest to the Fermi level for a given k. It evolves 
continously with temperature from zero temperature up 
to very high temperatures where there is no sharp peak in 
the one-particle local denisty of states (LDOS, Fig. 1(e)). 
The character of the dispersive excitations evolves conti- 
nously from strongly renormalized qps located near the 
Luttinger Fermi surface with Fermi crossings around the 
X point and on the F-M line for T <C Tsat (Fig. l(a-b)) 
to holes in the lower Hubbard band (located near the M 
point) (Fig. 1(d)) for T :^ Tgati as the spin degrees of free- 
dom gradually unbind from the charge, with increasing 
temperature. The qp velocity is nearly temperature in- 
dependent only below TpL and above Tgat- The mass en- 
hancement (1/^), decreases with increasing temperature, 
from a large value ~ 5 below TpL (Fig. 1(a)) to a value 
~ 1.5 c^ (1 — n/2)~^ at high temperature (Fig. 1(d)). 

We now turn to the transport properties and focus on 
the electric current induced by electric fields and thermal 
gradients. 



implementation of ref. [22] to solve the auxiliary impu- 
rity problem. We use Fade approximants to analytically 
continue the self energy. 

The one-electron spectral function is defined as 
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in terms of the bare band dispersion ek 



-1(C05(/C^ 



cos{ky)) and the self energy T^{uj). A\^{uj) at different 
temperatures are plotted in Fig. 1 (a-d). 

Several characteristics of the evolution of A\^{uj) with 
temperature are important. The solutions of the follow- 
ing equation. 



uj + i^{T) - ek - KS(cj, T) = 0, 



(3) 



faithfully reproduce the location of the peaks and how 
they evolve with temperature in A\^{uj) (Fig. 1(f)). We 
do not describe in the following the upper Hubbard band, 
at positive energies of order U. 

There are two distinct temperature regimes separated 
by a crossover scale Tgat — '^Tbr/S = 0.1, which also sets 
the scale for the saturation of qp scattering rate as will 
be explained later. Above Tgat^ say at T = 0.125, A\^{uj) 
has one peak, i.e., Eq. 3 has only one root for each k 
and displays a continuous dispersion over the whole Bril- 
lioun zone. Below Tsat^ Eq. 3 can have multiple roots. 
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a^/^ is charge/thermal conductivity matrix. Several 
quantities of interest are resistivity (p). Hall angle (6^7:/), 
Seebeck coefficient (5*) and Nernst coefficient (i^) [20]. 
They are representative measures of longitudi- 
nal/transverse and magneto/thermo-electric transport 
properties and can be expressed in terms of elements of 
conductivity matrices. 
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Within the DMFT treatment of the one-band Hubbard 
model current vertex corrections vanish and the trans- 
port properties can be interpreted directly in terms of 
one-electron spectral function [5, 8, 9], 
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with Of = or 1 for charge or thermal conductivity. We 
consider the hmit of weak magnetic field, hence the off- 
diagonal conductivities are proportional to B. ^ 
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are transport functions in 



terms of bare band dispersion ek and its derivatives. The 
derivatives are denoted by corresponding superscripts, 

To recast Eqs. 6 in terms of qp , we linearize Eq. 3 at 



and define Zk = (1 
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Then the low 



energy part of the one-electron Green's function can be 
approximated as 
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Thus Zk is the qp renormalization factor (or qp weight) 
and r^ = — Zk^5](cc;^) is the qp scattering rate. 

Then the integrals in Eqs. 6 can be performed analyt- 
ically and lead to 
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(r^) ^ is the qp lifetime. Transport functions 



are renormalized by Zk. ^k^ 
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This reformulation leads to a transparent interpreta- 
tion in terms of qps with temperature dependent disper- 
sion cj^. Eqs. 8 has a form similar to the solution of the 
kinetic equations from Boltzmann theory [23]. The es- 
sential difference from the Prange-Kadanoff treatement 
of the electron-phonon problem [2] is the strong temper- 
ature dependence of the qp dispersion brought in by Zk. 

First we validate the simplified description of trans- 
port, Eqs. 8 ("qp approx."), by benchmarking it against 
the result of the exact DMFT calculation, Eqs. 6("full 
calc"), for the resistivity. Hall angle, Seebeck and Nernst 
coefficients. The quantitative agreement between Eqs. 8 
and Eqs. 6 is evident, as shown in Fig. 2(a)-(d). The 
qp approximation faithfully reproduces the results of all 
transport quantities over the whole temperature range, 
extending to temperatures well above Tgat- 

Fig. 3(a) shows the qp scattering rate on Fermi sur- 
face, i.e., F^ with uol. =0 (For later use we also write 
r^ = (F^ )~^ as the qp hfetime and Zk^ as the renor- 
malization factor at Fermi surface). Tgat demarcates 
the non-monotonic temperature dependence of F^ . Be- 
low Tsati r^ increases and reaches maximum at Tgat- 
Above Tsat^ r^ decreases very slowly and eventually 
approaches to a value moderately smaller than the maxi- 
mum. This confirms that Tgat characterizes the crossover 
between two distinct scattering behaviors. The inset of 
Fig. 3(a) shows estimated values of {kpl*)'^ with kp an 
estimation of the average Fermi momentum by assuming 



a circular Fermi surface containing (1 — 5)/ 2 electrons 
and with /* = v^ r^ the qp mean free path, where 

'^k ~ \/('^k) ^^^^ (• • • ) averaging over the Fermi level. 
At low temperatures, {kpl*)'^ increases with tempera- 
ture, as expected in a good metal, and crosses over to 
a much slower increase, or saturated behavior around 
T5at/2. Above Tsat/"^, {kpl*)'^ ^ 0.5, and does not 
exceed the Mott-Ioffe-Regel (MIR) bound, which states 
that ikpl*)'^ < 1 in a metal. The qps behave as ex- 
pected in Boltzmann transport theory in the full tem- 
perature range, reaching the non-degenerate limit at 
T ^ Tsat' Notice that above TpL^ ^E(O) is not quadratic 
in temperature, only F^ = — Zk^^E(O) is quadratic. 

The anomalies in the transport properties are the re- 
sult of the strong temperature dependence of the renor- 
malized dispersion. This is best understood by means 
of a general Sommerfeld expansion of Eqs. 8, which is 
explained in the Supplementary Material and works well 
below Tsat 12. For this purpose we define <|)*^^/^^(e) = 

X^k ^k^^ ^^^(e — c^k) ^^^ ^^ energy dependent qp life- 
time r*(e) = r^ when e = cj^, with scattering rate 
F*(e) = (r*(e))-^ For |e| < T, $*^^/2/^(e) is expanded 
to the linear order in e. To keep the asymmetry in F*(e) 
which is important for the thermoelectric transport, we 
expand F*(e) to cubic order of e, and treat the linear and 
cubic order as corrections to the zeroth and quadratic 
terms, which are dominant in the Fermi liquid regime 
at low temperatures. The insets in Fig. 2 compare the 
estimation using this expansion (purple dots) and the re- 
sults of the full expressions (black dots). The agreement 
is evident and the expansion quantitatively captures the 
variation below Tsat 1 2. 

The inset of Fig. 2(a) shows the linearity of resistivity, 
a typical non- Fermi liquid behavior [24], up to Tgat/4 — 
0.025, as indicated by the linear fitting (blue dashed line). 
Surprisingly the qp scattering rate FJ^ has a quadratic 
temperature dependence also up to Tsat/^ (Fig. 3(b)). 
This is due to the strong temperature dependence of Zk^^ 
(Fig. 3(c)). In fact, the leading order in the expansion 
gives 



(Zk,$--(/i)r^J- 



(9) 



where <l>^^(/i) = Y^^ <l>^^5(/i-ek) with jl = /i-KE(O) and 
we have used <l>*^^(0) = Zk^<^^'^(/i). Zk^ c^ 0.1 + 12T 
for TpL <T < Tsat 1^1 leads to the quasilinear resistivity 
and also affects all other transport coeffecients in Eqs. 8. 
The temperature dependence of Zk^ becomes negligible 
only below the Fermi liquid temperature TpL — ^sat /15. 
$^^(/i) is very weakly temperature dependent as shown 
in the inset of Fig. 3(c). Above Tg^t, the resistivity is 
quasilinear in temperature with a slope smaller than that 
below Tsat 1^1 while the qp scattering rate is saturated. 
The expansion cannot be used at high temperatures, but 
the discrepancy between the scattering rate and the re- 
sistivity, can be traced again to the temperature depen- 
dence of the dispersion, which in this regime is due to 
the the shift of jl with temperature and not to changes 
in qp band structure. 
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FIG. 2: Transport properties, (a) Resistivity, (b) Hall angle, (c) Seebeck coefficient, (d) Nernst coefficient. Black dots ("full 
calc") are obtained using Eqs. 6. Red dots ("qp approx.") are obtained using Eqs. 8. Purple dots ("expansion") are obtained 
using Sommerfeld-like expansion which is detailed in Supplementary Material. 
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FIG. 3: (a) Quasiparticle scattering rate F^^. The inset shows the estimation of (kpl*) ^- (b) FJJ.^ as a function of T^ for 
T < Tsat/"^- (c) Quasiparticle renormalization factor Zi^^ for T < Tsat/'^- 



Similarly, the leading order in the Hall angle (Fig. 2(b)) 
is given by id^nOH/B c^ Z\^p.^y^{jl)r^ /^^^{jl) and indi- 
cates the sign change at T :^ Tgat/^ is due to the sign 
change in ^^^{jl), a consequence of the evolution of Fermi 
surface from a hole-like one to an electron-like one. For 
the Seebeck coefficent (Fig. 2(c)), the expansion leads to 
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The asymmetry in scattering rate competes with the 
asymmetry in the qp band structure, hence instead of 
sign change, 5* shows non-monotonic temperature depen- 
dence below Tsat/"^' 

The Nernst coefficient u (Fig. 2(d)) rises steeply below 
Tsat/^1 and provides a good probe of the temperature 
dependence of r^ . The leading orders in the expansion 
gives 
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(11) 
In the square lattice near hall-filling, the asymmetry in 

band structure dominates and leads to i^ ex r^ T ex 1/T. 
This rise is seen in many materials [20] before v drops 
linearly in T at very low temperature [37]. 

Further studies should be carried out to ascertain to 
which extent the DMFT description of transport ap- 
plies to real materials, but the strong similarities be- 



tween the experimental features revealed in the phe- 

nomenological picture in ref. [24] and our results are 

encouraging. AC transport measurements can be used 

to extract the temperature dependence of r^ . At 

low frequency, the optical conductivity is parametrized 
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as [25, 26] a{uj) = 

in AC Hall effect [27], B tan Oh{uj) = 
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and tI 



T^ /2. Frequency dependent thermoelectric measure- 
ments would give additional information on the asym- 
metry of the qp dispersion and scattering rate. 

The extension from model Hamiltonians to the 
LDA+DMFT framework is straightforward. It can be 
used to separate the temperature dependence of trans- 
port coefficients arising from the temperature depen- 
dence of the qp band and that of the scattering rate, in 
materials such as the ruthenates [28], the vanadates [29] 
and the nickelates [14, 30, 31] for which the LDA+DMFT 
description is known to provide an accurate zeroth order 
picture of numerous properties [14]. Recent experiments 
on cuprates [32] have revealed evidence for temperature 
dependence of cj*^^ and a T^ scattering rate over a broad 
range of temperatures. These materials require cluster 
DMFT studies to describe their momentum space differ- 
entiation. Still, it is tempting to interpret the transport 



properties in terms of qps to provide an effective de- 
scription of the transport. Indeed the qp scattering rate 
computed in the t-J model in ref. [33], exhibits the sat- 
uration behavior described in this work and it would be 
interesting to re-analyze the results in terms of the qps 
of the hidden Fermi liquid. Our findings are related to 
two earlier theoretical proposals. Anderson introduced 
the idea of a hidden Fermi liquid [34, 35], requiring Zk^ 
strictly vanishing at T = in the normal state. Alterna- 



tively, our results could be cast into the framework of the 
extremely correlated Fermi liquid [36] by the temperature 
dependence of the caparison function. 
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In this Supplementary Material, we consider the low 
temperature behavior of the integral 



~ J T \t) 
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F{e) 



4cosh^{^)G{e)' 
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In evaluates the conductivities when F{e) and G(e) are 
chosen properly. For F{e) = <l>*^^(e) and G(e) = r*(e) = 
(r*(e))~^, /o gives the longitudinal conductivity cr^^ and 
/i gives the thermal counterpart a^.^. For F{e) = <I>*^^(e) 
and G(e) = (r*(e))^, /q gives transverse conductivity ay^ 
and /i gives ciy^. 

At low temperature, <l>*^^(e) and <l>*^^(e) are regular 
functions that vary smoothly for |e| < T. Thus for the 
purpose of evaluating the intergral, they can be expanded 



as, 






$*2/^(0). 



• <^*^^ (0)e, 
<l>*^^'(0)e. 



(2) 
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where <l>*^^ (0) = d<l>*^^(0)/de and <l>*^^ (0) = 
d<l>*^^(0)/de. 

The quasiparticle scattering rate r*(e) in a Fermi liq- 
uid is dominated by the particle-hole symmetric part at 
low temperature, and has the quadratic form 



r*(e)^5o + 52e' 



(4) 



In a clean system, go ex T^ and consequently, r*(e)~^ be- 
comes a Lorentzian function with vanishing width when 
T ^ 0. Hence r*(e)~^ is asymptotically singular for 
|e| < T and the standard Sommerfeld treatment is no 
longer appropriate [1]. Furthermore, the particle- hole 
asymmetric part in r*(e) is in particular important for 
thermoelectric transport. The asymmetric part can be 
approximated as the linear and cubic order in e as cor- 
rections to the quadratic scattering rate [2]. Then we 
write 



r*(e) c^go^g2e^^gie^g3e^. 



(5) 



This also leads to the corrections to the Lorentzian form 
of r*(e)~^, which are expressed as 
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By changing the variable e = xT, we define the follow- 
ing integral. 
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^ 4cosh\x/2) {go/T^ ^ g2X^)^' 
and the longitudinal conductivity is estimated by 
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-g^T^*^*^^' {0)Eiy 



(9) 



The first term is determined by the leading quadratic 
term in Eq. 6. The second and third terms are deter- 
mined by the linear and cubic corrections, i.e., propor- 
tional to gi and gs respectively. Similarly, for other con- 
ductivities, we have 
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In each expansion, the second and third term will be 
simply denoted as the gi and gs term. 

In general, go, ^i, ^2 and gs vary with temperature, 
thus E^^ also depends on T in a non-trivial way. But for 
a Fermi liquid without impurity scattering, when T ^ 0, 
(^0 + ^'2^^) ^ (tt^T^ + e^), hence E^ is a constant inde- 
pendent of T. Besides, ^1 ex T^ as T ^ [2]. We obtain 
the low-temperature leading terms of conductivities, a^^, 
and ay^ are dominated respectively by the T~^ and T~^ 
divergence in the first term. The prefactors are deter- 
mined by the qp band structure (<l>*^^(0) and <I>*^^(0)). 
For the thermoelectric transport, all three terms lead to 
the same power law (T~^ for a^^ and T~^ for Cy^)^ and 
the prefactors are determined by the asymmetry of qp 
band structure (<l>*^^ (0) and <l>*^^ (0)) and scattering 
rate {gi and ^3). 
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FIG. 1: Conductivities computed using Eq. 1 ("qp approx." 
and using the expansion of Eq. 9 — 12. 



The Sommerfeld expansion is equivalent to neglecting 
the g2e^ term in the expansion of scattering rate, Eq. 6 
and Eq. 7. Above TpL^ 92 is small, and the Sommerfeld 
expansion can give quantitative estimation that correctly 
captures the variation of conductivities with tempera- 
ture. The Mott relations for thermoelectric transport, 
i.e., 



are based on the Sommerfeld expansion [3, 4], thus they 
still provide convenient rules for estimating the Seebeck 
and Nernst coefficient when Tfl < T < Tsat/"^- Be- 
low TpL, the g2e^ term will lead to different prefactors 
(Eq. 8), but the Fermi liquid power law at low tempera- 
ture is not changed. 



We compute conductivities for T < Tgat/'^ using Eq. 1 
("qp approx.") and using the expansion of Eq. 9 — 
Eq. 12. The results are compared in Fig. 1. The agree- 
ment between the two sets of data is evident. Therefore 
Eq. 9 — Eq. 12 provide a convenient way to decompose 
the transport properties. In our calculation, a^^ is dom- 
inated by the first term of the expansion. Hence the 
resistivity is determined by <^*^^(0) and the quasiparti- 
cle scattering rate at the Fermi level r*(0) (equivalent to 
Fk^ defined in the main text). Similar observation can 
be drawn for a^^. In a^^,, the g^ term is negligible but the 
first and the gi term have opposite signs and both change 
sign below Tsat/"^- This is due to the different temper- 
ature dependence of the asymmetry of the quasiparticle 
band structure and scattering rate and they together de- 
termine <j^^ and thermopower. (jy^ is dominated by the 

first term because <l>*^^ (0) is close to its singularity at 
half-filling. 
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